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Recent results of the KASCADE air shower experiment are presented in order to shed some light on the 
astrophysics of cosmic rays in the region of the knee in the energy spectrum. The results include investigations of 
high-energy interactions in the atmosphere, the analysis of the arrival directions of cosmic rays, the determination 
of the mean logarithmic mass, and the unfolding of energy spectra for elemental groups. 



1. Introduction 

The origin of high-energy cosmic rays is among 
the most interesting questions in astrophysics. 
The origin of a structure in the all-particle en- 
ergy spectrum around 4 PeV, the so-called knee, 
is generally believed to be a corner stone in the 
understanding of the astrophysics of high-energy 
cosmic rays. The knee is proposed to be caused 
by the maximum energy reached in cosmic-ray 
accelerators or due to leakage of particles from 
the Galaxy. Hence, an understanding of the ori- 
gin of the knee reveals hints on the acceleration 
and propagation of cosmic rays. Experimental 
access to the understanding of the sources, accel- 
eration and propagation mechanisms is provided 
by detailed investigation of the arrival directions, 
energy spectra, and mass composition of the ul- 
trarelativistic particles. 

While at energies below 1 PeV cosmic rays can 
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be measured directly at the top of the atmo- 
sphere, the strongly decreasing flux as function 
of energy requires large acceptances and exposure 
times for higher energies. Presently they can be 
realized in ground based facilities only. There, 
the secondary products, generated by interactions 
of high-energy cosmic-ray particles in the atmo- 
sphere, the extensive air showers, are registered. 
It turns out that a correct description of the high- 
energy interactions in the atmosphere is crucial 
for a precise astrophysical interpretation of air 
shower measurements. 

To investigate cosmic rays from several lO''^ eV 
up to beyond 10'^ eV the air shower experiment 
KASCADE ("Karlsruhe Shower Core and Array 
DEtector") [1] is operated since 1996. The ex- 
periment detects the three main components of 
air showers simultaneously. A 200 x 200 m^ scin- 
tillator array measures the electromagnetic and 
muonic components. The 320 m^ central detec- 
tor system combines a large hadron calorimeter 
with several muon detection systems. In addition, 
high-energy muons are measured by an 128 m^ 
underground muon tracking detector. 
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2. High-energy interactions in the atmo- 
sphere 

A correct understanding of high-energy inter- 
actions in the atmosphere is indispensable for a 
good astrophysical interpretation of air shower 
data. The electromagnetic part of the showers 
is well understood and described by QED. For 
the air shower dcvclopmcuit the luiderstanding of 
multi-particle production in liadronic interactions 
with a small momentum transfer is essential. Due 
to the energy dependence of the coupling constant 
ag and the resulting large values for soft interac- 
tions, the latter cannot be calculated within QCD 
using perturbation theory. Instead, phenomeno- 
logical approaches have been introduced in differ- 
ent models. 

For the numerical simulation of the develop- 
ment of air showers the program CORSIKA [2] 
is widely used. It offers the possibility to use 
different models to describe low and high-energy 
hadronic interactions. A principal objective of 
the KASCADE experiment is to investigate the 
air shower development in detail and test the va- 
lidity of the models included in simulation codes 
such as CORSIKA, using as much information as 
possible from the simultaneous measurement of 
the electromagnetic, muonic and hadronic com- 
ponents. 

With these investigations already several prob- 
lems in existing codes could be pointed out and 
some interaction models (or particular versions of 
them) could be shown to be incompatible with the 
measured data [3-6]. As an example of present 
activities. Fig. 1 shows the hadronic energy sum 
as function of the number of electrons. Presented 
are measured values compared to predictions of 
three different interaction models for showers in- 
duced by primary protons or iron nuclei. For a 
presentation of the data as function of the num- 
ber of electrons one expects an enrichment of light 
primaries within the particular intervals, hence, 
the data should approach the values for the pro- 
ton component. One recognizes that the data are 
compatible with the predictions of QGSJET and 
SIBYLL, while on the other hand NEXUS 2 pre- 
dicts too less hadronic energy in most of the elec- 
tron number range. From such distributions one 
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Figure 1. Hadronic energy sum versus number of 
electrons. 

can conclude that the present version of NEXUS 
is not compatible with the data. More detailed 
investigations of the models QGSJET 01 and 
SIBYLL 2.1 are presently in progress. 

3. Arrival directions 

The investigation of the arrival directions of 
cosmic rays improves the understanding of the 
propagation of the particles through the Galaxy 
and their sources. Model calculations of the dif- 
fusion process in the Galactic magnetic field indi- 
cate that there could be an anisotropy on a scale 
of 10~^ to 10~^ depending on particle rigidity 
as well as strength and structure of the Galac- 
tic magnetic field [7]. Diffusion models relate a 
rigidity dependent leakage of particles from the 
Galaxy to the steepening (or knee) in the all- 
particle energy spectrum around 4 PeV. Thus, 
anisotropy measurements can provide substantial 
information on the origin of the knee. 

In KASCADE investigations [8,9] attention has 
been drawn to a search for point sources and 
large-scale anisotropy. Of special interest is the 
search for potential gamma-ray induced showers. 
Since photons are not deflected in the Galac- 
tic magnetic field, they arc suitable for a direct 
search for the sources of high-energy particles. 
Experimentally such an investigation is realized 
by the selection of muon poor showers. 

A search for point sources was performed for 
primary photon candidates as well as for all 
(charged) particles. The search covers the whole 
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Figure 2. Rayleigh amplitudes A as function of pri- 
mary energy. The KASCADE upper limit (bold line) 
is compared to results from the literature [10-13]. 
Model predictions [7] for the total anisotropy as well 
as for the light and heavy component are also shown 
(thin lines). 

sky, visible by KASCADE (declination 15° to 
80°). Special attention was drawn to the Galac- 
tic plane and known gamma-ray sources in the 
TeV region. The search was complemented by an 
analysis of the most energetic showers registered 
and an investigation of the most photon-like pri- 
mary particles. None of the searches reveals an 
indication for a significant excess of the flux. The 
distributions of the significance values over the 
whole sky follow the expectations for an isotropic 
cosmic-ray flux. 

The results for the analysis of the large-scale 
anisotropy are illustrated in Fig. 2. There, the 
Rayleigh amplitude is plotted versus the primary 
energy. The upper limit derived from the KAS- 
CADE data is compatible with other measure- 
ments from the literature and also with theoreti- 
cal calculations [7] for different elemental groups. 

4. Energy spectra for elemental groups 

The most direct experimental access to the as- 
trophysics of cosmic rays is provided by measure- 
ments of their energy spectrum and mass com- 
position. Investigations of the KASCADE ex- 
periment reveal that the knee in the all-particle 
spectrum is caused by a turn-off of the light com- 



ponent (protons and helium nuclei) [14]. This 
is implies an increase of the mean logarithmic 
mass {{In A)) as function of energy. While it is 
beyond doubt that the data reveal such an in- 
crease, the absolute values of (In^) depend on the 
observables investigated and on the interaction 
models used in the simulations to interpret the 
data [15-17]. The differences amount to about 
A (In A) « 1. This indicates that further and 
more detailed investigations of high-energy inter- 
actions in the atmosphere are necessary for an 
unambiguous astrophysical interpretation of the 
observed data. 

Presently, the most promising approach is the 
unfolding of energy spectra for individual elemen- 
tal groups from the data of the electromagnetic 
and muonic component [18-21]. Systematic stud- 
ies are performed with different hadronic interac- 
tion models available in CORSIKA, for both, low 
and high-energy (> 80 GeV) interactions. Dif- 
ferent unfolding methods are applied, in order to 
investigate the systematic effects introduced by 
the individual methods. The different methods 
result in similar flux values if the same code is 
used for the air shower simulation. But significant 
differences occur for different interaction models. 
This is illustrated in Fig. 3, where recent results 
of an analysis applying the Gold algorithm are 
depicted. 

The figure also shows the flux of primary pro- 
tons as obtained in an analysis of unaccompanied 
hadrons [22]. The flux is compatible with the 
proton flux as obtained from the unfolding pro- 
cedure. For comparison, also the results of direct 
measurements at lower energies at the top of the 
atmosphere [23] are given in the flgure. 

The analyses indicate that, at present, the 
understanding of primary cosmic rays is lim- 
ited by the insufficient knowledge of high-energy 
hadronic interactions in the atmosphere and not 
by too low statistics of the measurements or the 
systematic uncertainties introduced by different 
reconstruction methods. 

5. Conclusion and outlook 

The present status of the KASCADE exper- 
iment to investigate high-energy cosmic rays in 
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Figure 3. Energy spectra for groups of elements derived from the data using CORSIKA with the hadronic 
interaction models QGSJET (left) and SIBYLL (right). The lines indicate extrapolations of direct measurements 

[23]. 



the knee region has been discussed. Continuing 
efforts are taken in order to improve the under- 
standing of high-energy interactions in the atmo- 
sphere. A search for point sources of charged par- 
ticles as well as gamma rays did not reveal any 
significant excess. Upper limits for the large scale 
anisotropy have been derived. The observations 
reveal an increase of the mean logarithmic mass 
as function of energy in the knee region. Energy 
spectra for groups of elements have been derived. 
To extend the investigations to higher energies up 
to 5 X 10^^ eV data taking with an enlarged array 
[24] has started in 2003. 
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